Exact Solution of Ising Model on a Small- World Network 
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We present an exact solution of a one-dimensional Ising chain with both nearest neighbor and 
random long-range interactions. Not surprisingly, the solution confirms the mean field character 
of the transition. This solution also predicts the finite-size scaling that we observe in numerical 
simulations. 
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I. INTRODUCTION 

Physical interactions are usually of limited range. Nev- 
ertheless, there is a long history of study of the effect 
of infinite range interactions in basic models of statisti- 
cal physics P, |3| • More recently, the work of Watts and 
Strogatz Q on small world networks brought renewed at- 
tention to this question. Watts and Strogatz considered 
a quantitative model that interpolates between a regular 
lattice and a random one. They showed that the addition 
of a small fraction of random long range links can dra- 
matically change the connectivity properties of the lat- 
tice. In particular, in a lattice of N nodes the averaged 
chemical distance, that is the averaged minimum num- 
ber of links between any two nodes, is of order 0{hiN), 
rather than 0{N), as in a regular lattice, regardless of 
the number pN of additional Hnks (p > 0), provided TV 
is large enough. Clearly, this change should be reflected 
in the phenomenology of any physical models deflned on 
such lattices. Researchers were quick to grasp this op- 
portunity and many physical models and processes have 
been considered in small-world networks, like cellular au- 
tomata Qjdiffusion 0, neural networks @, the spread 
of disease 0, and many others 

The Ising model, the simplest paradigm of order- 
disorder transitions, has been studied on small-world 
networks in one dimension (ID) both numerically and 
analytically H EHH El, and in 2D and 3D numeri- 
cally Some partial analytic results, based on the 
use of the replica trick, were presented by Gitterman 
[m and Barratt and Weigt E3|- They are, however, 
contradictory (Gitterman predicts a ferromagnetic phase 
only for p > 1/2) and, at any rate, these results do not 
constitute a full solution of the thermodynamics of the 
model. There is a consensus that the ferromagnetic Ising 
model has a mean-field transition for any finite value of 
p even in ID, even though Hastings predicts the occur- 
rence of anomalous scaling of the mean-field amplitudes 
with p — > E3| • Dorogovtsev et al E3| presented a very 
general discussion, for networks that are locally tree-like, 
based on the solution of the Ising model on a Bethe lat- 
tice 

In this work we present an explicit solution of the ID 



Ising model with additional random long-range bonds. 
We are able to calculate not only the transition tem- 
perature, but also the complete thermodynamics at all 
temperatures as a function of p and of the strength of 
long range interactions / (though in this paper we only 
present the results for the free energy and specific heat, 
in zero field) . The form of finite size scaling corrections 
to the thermodynamics can also be derived from our so- 
lution. 

The following Sec. lUlpresents a description of our spe- 
cific model and the basic algorithm for its exact solution. 
The combinatorial treatment of its constituents is given 
is Sec. mil In Sec. IIVI we analyze the resulting ther- 
modynamics and conclude on the mean-field character 
of the ordering transition. Sec. shows how the finite 
size scaHng properties can be obtained within the same 
approach. Finally, the relation between our results and 
those of Dorogovtsev et al is discussed in Sec. IVII 



II. THE MODEL 

Watts and Strogatz originally considered a model in 
which the bonds of a regular lattice are rewired at ran- 
dom with a probability p. It is widely believed that the 
modified model, in which random long range bonds, or 
shortcuts, are added to the regular lattice |3, is essen- 
tially equivalent. 

The actual model consists of a chain of N Ising spins 
Ui, with nearest neighbor interactions J (chain bonds) 
and shortcut interactions / (long range bonds), both J 
and / being positive, so that the Hamiltonian reads 
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with (Ti = ±1 and periodic boundary conditions: ctq = 
ctat. The set S contains Nb = pN shortcut pairs of spins, 
and the last term accounts for the effect of external mag- 
netic field. We present in detail the solution of a version 
of this model in which the sites connected by shortcuts 
are equally spaced in the regular lattice, a distance l/2p 
apart, but the shortcuts are randomly arranged among 
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these spins. The solution of theoriginal model, in which 
the 2pN sites are randomly distributed along the lattice, 
turns out to be essentially the same, and some results 
for this latter case are also presented in the Appendix. 
In either case, the bond selection is such that no spin is 
linked to more than one shortcut. 

We consider a transformation from site variables to 
bond variables, which we formulate for a general Ising 
model in zero magnetic field, 

n = - ^JijUiaj, (2) 

where Jij can be chosen arbitrary (though in what fol- 
lows they are only 0, J or /). The partition function for 
temperature T = /?~^ is given by 

Z = Tr{^}exp(/3^ Jy cr,;(Tj) (3) 

= Tr{^} Jl CXp(/?Jy CTjCrj). 
(ij) 

We can use the known identity 

exp(pjijai<7j) — coshpjij (1 + ai<7j taiih fJJij) (4) 

and then represent the latter factor as a sum in a discrete 
"bond variable" bij = 0, 1: 

1 + (TiCTj tanh/3Jy = {cTicrj ta,nh (3 Jij)^'\ (5) 

6ij=0,l 

to write the partition function as 

Z = J]cosh/3Jy (6) 

\ih3) J 

xTr{,} n (^^'^J tanh/3Jy )''- . (7) 
{b} 

When we trace over any spin variable ai, in a product 
with fixed configuration {b} of bond variables bij, we get 
zero if the ith spin multipHcity in this product, bij, 
is odd and a factor of 2 if it is even. Therefore we can 
trace over the spin variables to obtain 

Z = 2^ I J] cosh/3Jy I (8) 
X E n (tanh/JJ,,)"- , (9) 

where the sum over {b} is restricted to configurations 
with only even J^j bij for all i. 

As an example, we can derive from Eq. ^ the partition 
function of the Ising chain {Jij = J, for nearest neigh- 
bors). Since each spin there has only two bonds, this 
implies that either bij = for all the bonds or bij = 1 



for all the bonds, and these two alternatives contribute 
in the partition function 

Zcha.n = (2 cosh f3jf (l + tanh^ j) (10) 

as expected. 

The restriction of the sum in Eq. Q to configurations 
such that J2j ^ij is even, allows a reduction of the number 
of required variables. If the site i has no shorcut bonds, 
then this condition reduces to + being an even 
number, that is, bi-i,i = (since bij = 0,1). On 

the other hand, if i is a shorcut site, with an associated 
additional bond variable, the sum in Eq. ^ gets only 
nonvanishing contributions when = if the 

additional bond variable is 0, or when = 1 — 

if the additional bond variable is 1. In any case, the 
knowledge of one chain bond variable (say 6o,i) and the 
shorcut variables is enough to determine all the terms 
that give a nonzero contribution to Eq. Hence, the 
partition function can be written as an unrestricted sum 
over 6o,i and iVf, shorcut bond variables which, in order 
to simplify notation, we denote simply by 6o, 5i, . . . , bN^: 

Z = Z^,^,^cf Yl 'J^'^'f^ (11) 

where cj = cosh/3/, tj = tanh/3J, tj = tanh/3/, 
M — X)r=i is the number of shortcut bonds with ^ = 1, 
and L the number of chain bonds with 6 = 1. We calcu- 
late it as follows: given a configuration of iVf, shortcuts 
and a choice of M from them with bi = 1, the chain gets 
divided into 2M segments between consecutive spins con- 
nected to one of M bonds. These segments have lengths 
h, ■ ■ ■ , hhi- When bo = 0, the value L[b] is the sum of 
even lengths, l2 + li + ■ ■ ■ + l-iM- Otherwise, it is the sum 
of the odd ones, h + I3 + . . . + I2M-1. 

Our solution is based on the observation that Eq. I|ll|) 
permits explicit separation of two extensive thermody- 
namical variables L[b] and M[b] (both being 0{N)) from 
N(, microscopic variables bi. The latter only define a tem- 
perature independent prefactor fl{M, L) at the product 
tjtf'' with given L and M. A similar situation (but with 
a single variable tj) was already addressed by Saul and 
Kardar for the problem of Ising spin glasses I^" 
stead of numerical procedure for the prefactor employed 
in Ref. are able here to calculate analytically the 

sum over all the bond configurations in Eq. ifTTll . Denot- 
ing fl{M, L) = exp[S'(Af, L)] the number of choices of M 
segments with a fixed value of L, we get 

bo M,L 

where ^,7 = l/ln(l/tj) is the Ising chain correlation ra- 
dius and ^/ = l/ln(l/t/). 

The crucial property of the exponential function un- 
der sum in Eq. lfT2|l is that it has, as will be seen below, 
an extremely sharp maximum at some macroscopically 
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great values of L and M. This reduces the calculation of 
its contribution to the free energy, = — TlnZ, to the 
logarithm of its maximum. The latter turns out to be 
0{N) and thus a self-averaging quantity, that is coincid- 
ing, for (almost) any random realization of the disordered 
system, with its average value (T^ . 

Passing to the intensive thermodynamical variables I = 
L/N, n = M/Nb, and to the intensive function s{l,n) = 
S{L, M)/N, we rewrite Eq. (O as 

Z = 2Z,,,„„c?^ ^^P[^ - ^Z^'' - P"/^/)]' (13) 

and in the thermodynamic limit N 00 the sum con- 
verts into integral which can be done by steepest de- 
scent. The values 60 = and 1 give identical contri- 
butions, hence the factor 2. The free energy per spin is 

/ = J^/N = -{T/N) InZ = -{T/N) InZ^Za = h + fa, 
where 

/o = -T[ln(2 cosh (3J) + p In cosh [31] (14) 

and the "anomalous" term fa (if exists, see below) is given 
by simple minimization with respect to I and n of the 
following function 

f{l,n) = T[-s{l,n) + l/tj+pn/^i]. (15) 

The task that remains is to calculate the configurational 
entropy s{l,n) = \nn{L,M)/N (for L,M = 0{N)). 



III. COMBINATORICS OF BONDS 

The calculation of ri(i, M) can be formulated in the 
following way. 

We are given a chain of TV sites {0, 1, 2, . . . , TV — 1} 
with periodic boundary conditions. The shortcut sites 
are evenly spaced forming a regular lattice with coordi- 
nates {d, 2d, . . . , TV}, where d = l/2p. A number pN of 
shortcuts connect pN pairs, randomly chosen from these 
sites. If we choose AI of these bonds (those for which 
6=1) from the total of pTV, the corresponding shortcut 
sites (the filled dots in Fig. iQJ) will have coordinates, 
in increasing order {ri, r2, . . . , r2A/} ^ {d, 2d, ... , TV}. 
These site coordinates will divide the lattice into 2TW 
segments of lengths dli, . . . , dl2M where 

k = '-^^^^ * ^ 1: (16) 

TV - r2M + n 



and — N/d. Then ri{L,M) is the number of 

possible choices of the M bonds such that 

l2 + h + k + ■■■ + I2M = L/d. (17) 




Figure 1: Q{L,AI) is the number of possible choices of M 
bonds, terminating in filled dots, from pN fixed bonds, which 
have a given value of the sum h + U + ■ ■ ■ + I2M = L. 

We are asking in how many ways one can di- 
vide the interval [Q,N/d\ into 2M + 1 integer lengths, 
{?!,..., Z2M+1}, where only I2M+1 may be zero, so that 
even lengths are summed to 

l2 + li + ■■■ + I2M = L/d. (18) 
Then, clearly, the odd lengths should sum to 

TV - L 

^1 + ^3 + • ■ • + I2M+1 = ^ ■ (19) 

Therefore Vl{L, M) is just the number of ways of dividing 
the interval [0,L/(i] into M segments of integer length, 
times the number of ways of dividing [0, (TV ~ L)/d] into 
M + \ integer segments (the last of them possible to be 
zero): 

Cir''"ciit-' (20) 

(where = n\/[m\[n — m)!], a binomial coefficient). 

Note, however, that our choice of 2M coordinates 
{ri, . . . , r2A/} is constrained by the the fact that the orig- 
inal problem has bonds connecting pairs of sites qk,qm- 
Therefore, if qu G {ri, . . . , r2M} so must qm. To take this 
fact into account, we multiply the previous factor by the 
normalizing factor, 

Cl7/C?M, (21) 

which refers to the self-averaging property in the ther- 
modynamic Hmit: that (almost) all possible pairings be- 
tween the sites {d, 2d, . . . , TV} give the same ^}{L,M). 
Considering Eq. Illl|l . we express 

niL,M) = Clf''-''^Clf_-'C^,f/C',lf (22) 

At this point it is important to specify the variation 
range for the variables L and M. Clearly, < TVf < pTV. 
On the other hand, since each ri — Vi-i measures at least 
d = l/2p, we must have 2pL > M and pN - 2pL > M, 
i.e. the before defined intensive variables n — M/pN and 
I = L/N should belong to the triangle 

n<2l, n<l-2l < n < 1, (23) 

shown in Fig. (j^J. 
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Figure 2: Trajectories of the maxima of the exponent in 
En.llHfl in the domain {l,n}, as temperature varies from 
T = (central point), to (lateral vertices, bo = or &o = 1) 
at particular choices of parameters: 1) p — 0.01, / = 1; 2) 
p = 1/2, 7 = 1; 3) p = 0.1, / = 0.1; 4) p = 1/2, I = 0.1. For 
each temperature there are two equal contributions to free 
energy resulting from bo = 1,0. 




Figure 3: Histogram of occurrences of I values for several 
fixed values of n, obtained in multicanonical simulations with 
the entropy of Eq. \2(\^ as sampling distribution for a sample 
with p = 1/4 and iV = 512. 



For our purposes, all that will be required is the leading 
0{N) term in hiVl{L, M). Using Stirling's formula 



we arrive at 



with 



sil, n) 



In r2(L,M) = Ns{l,n) 



p [(2 ~ 20 In (2 - 20 + 2nn(2/) 
+ (1 — n) ln(l — n) — nlnn 
-{2l-n) ln{2l-n) 
- (2 - 2Z - n) In (2 - 2^ - n)] , 



(24) 



(25) 



(26) 



(dropping the terms O(lnA^)). 

In order to confirm the basic assumption of self- 
averaging, involved in the calculations leading to 
Eq. Il2(i|l . we performed multicanonical Monte Carlo sim- 
ulations on the variables I and n. For the sampHng prob- 
ability distribution of a state i (characterized by a set of 
values for the variables {b}) we used pi = l/Q{Li, Mi). 
The Monte Carlo move is given by: 

• From an initial state i, propose a new state j, dif- 
fering from i by the value of a randomly chosen 
bond (with probability l/Nj,). 

• Accept the new state with the usual probability 



miii{l,n{L,,M,)/n{L,,M,)). 



This procedure ensures (by detailed balance and ergod- 
icity of the algorithm) that the asymptotical sampling 
probability is pi . The asymptotical histogram of frequen- 
cies in an (/, n) point should, then, be given by. 



H{l,n) cx Q.a{L,M) 



1 



n{L,M) 



where fla{L,M) is the actual number of states, with 
given {L, M), of the sample that we are using in the sim- 
ulation. If 17a(L,M) = n{L,AI), the histogram should 
be fiat in all the range of {l,n). 

The corresponding histograms (shown in (Fig. Q, for 
a single disorder reaHzation) are indeed quite fiat , con- 
firming that the entropies we calculated are apparently 
exact in the thermodynamic limit and very accurate 
for the moderate tested sizes. Recall that fl{L, M) (x 
exp(iVs(Z, n)) and, for large iV, even small deviations 
of In ria(L,M) from S{l,n), would result in simulations 
witch do not cover the entire spectrum, as ours do. 



IV. THERMODYNAMIC BEHAVIOUR 

The necessary conditions for that the exponent in 
Eq. ifn^ is maximum: df{l,n)/dn ~ df{l,n)/dl ~ 0, 
define from Eq. Il26p the following equations for the co- 
ordinates n* and /* of this maximum 



n*(l — n*) 



(2 - 21* - n*){2l* - n*) 
(21* - n*){l - r) 



ti 



(27) 



l*{2~2l* ~n*) 
(d = l/2p). These equations are easily solved: 



I* 



(1 - tjy - ititj ' 



r 



1 



ti - f^j 
i + ti 



(28) 
(29) 



and, for tj, tj varying with temperature, these are para- 
metric equations for a certain trajectory of the maximum 
of the exponent in the n, I plane (within the admissible 
range, Eq. ^21^ ). At zero temperature, when tj = tj = 1, 
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it locates at /* = 1/2, n* = 1/2 (the central point in 
Fig. I2II , corresponding to the maximum configurational 
entropy s{l,n). As T increases, the maximum moves to- 
wards the origin n = I = (or to n = 0, ^ = 1), as shown 
in Fig. (jSJ. The trajectory reaches the very origin at a 
finite temperature T = Tc, given by the condition 



4(1 + 211) = 1. 



(30) 



It is important to notice that below this critical temper- 
ature, when the reduced temperature t = T/Tc—1 is non 
zero (t < 0), no matter how small \t\ is, the coordinates of 
maximum (Z*,n*) Eqs. Ij28[ I29|l . are also finite: l*,n* ^ 
\t\. Hence the corresponding numbers L* ~ l*N and 
M* = pn*N are macroscopic: L*,M* = 0{N). On the 
other hand, the width of the maximum, estimated from 
the second derivatives 9^s/9Z^|j, ^, ~ 9^s/9n^|j, ^, '-^ 

l/l*, is 0{^/l*/N) and tends to zero in the thermody- 
namic limit, that is the probability distribution in macro- 
scopic variables I, n tends to a (5-function with an ampli- 
tude which can only depend on macroscopic parameters 
of long range links (p and /) but not on their specific 
realization. This justifies the above used assumption of 
self-averaging for the probability distribution and relates 
it to the known self-averaging property for observable 
values (1^ . 

Above Tc, the maximum of the exponent in Eq. Ijl3|l 
goes away from the physical region, Eq. lf23jl . while its 
highest value in this region is zero, attained at the ori- 
gin. Hence the last factor of the Eq. (jl3ll turns 0(1), 
and it gives no contribution to the free energy in the 
thermodynamic limit. Therefore the free energy above 
Tc is simply /o, Eq. ifnt . but it gets an extra term below 
Tc'- fo + fa, where fa is given by Eq. Ijl5|l with I and n 
given by Eqs. (l28)l and ll29ll . These equations combined 
provide a complete description of the thermodynamics of 
the model at all temperatures in zero field. 

The dependence of transition temperature Tc, defined 
from Eq. Ij30p . on the system parameters p, J, and / can 
be easily analyzed in characteristic limits. 

If shortcut bonds are much stronger than chain bonds, 
then for any finite p and / ^ 00, critical temperature 
tends to a finite value: Tc J/(arctanh9~P). Since 
shortcut spin pairs in this Hmit should be considered 
as single spins, the above value defines also the transi- 
tion temperature for a random graph with connectivity 
4 made of Ising chains of length d = l/2p. In the limit 
of small concentration, p 1, this critical temperature 
Tc turns small compared to J, the energy scale for Ising 
chain: 



Tc 



2J 



ln[l/(pln3)]' 



(31) 



and such /-independent behavior holds as well for mod- 
erate shortcut strength (unless / is too small: / <C 
J/l ln(pln3)|). At last, in the Hmit where the shortcuts 
are much weaker than the chain bonds, pi ^ J, we have 




^ HI/ J) ° 



Figure 4: Behavior of the critical temperature (in units of J) 
as a function of I for several values of the concentration p of 
shortcuts 




Figure 5: Specific heat obtained by Monte Carlo simulations 
for 4 sizes of samples, with p = 1/4 and 7 = 1 in comparison 
with the theoretical curve {00). Inset shows the collapse of 
the 4 curves when plotted in function of tN^^^ . 



within logarithmic accuracy 



2J 



ln[Tc/(2p/)] ln{J/[p/ln(J/p/)]}- 



(32) 



The above relations define the system phase diagram 
in p, I, T variables, as shown in Fig. 

It is of interest to compare these formulas to the fi- 
nite critical temperatures, resulting from breaking down 
the Mermin- Wagner theorem for a 2D Heisenberg mag- 
net in presence of (small) anisotropy AJ <^ J: Tc 
Jz/ln(Tc/AJ) {z the coordination number). They can 
be also referred to the percolation threshold Pc for the 
one-state limit of Potts model ^^: Tc = 2J/ln(l —pc)- 

The non-trivial thermodynamics follows from the ob- 
servation that close to the critical point, < — t ^ 1, 
both variables I and n are 0(1^1), while the trajectory 
n{l) reaches the origin with asymptotic slope dn/dl 
4/(l/t/ + 2) < 4/3, that is, always within the triangle, 
Eq. Using Eqs. (j^Hl and the additional term 
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fa in the free energy can be simplified to: 

(2-2/*-n*)2 



fa = Tp In 



[1 - n*) {2 - 2l*y 



(33) 



and its leading terms in the critical region are clearly 
of order 0{t'^). Hence the specific heat per spin, c.„ = 
—Td'^f/dT'^, has a finite jump AC at the critical point. 
As an example we present in Fig. © the calculated exact 
specific heat for p — 1/4 and I ~ J, together with the re- 
sults of simulations on samples of various sizes. This be- 
havior permits to classify the considered transition as sec- 
ond order mean-field-like, with critical exponent a = 0. 
In the case of p ^ 1 with Tc given by Eq. (13 111 , the jump 
is proportional to p^ln^(pln3), and is also independent 
of the shortcut bond strength /. When the shortcuts 
are much weaker that the chain bonds, pi ^ J, with 
Tc by Eq. Il32p . we find the specific heat discontinuity 
proportional to p{I / \n'{J/pI). The above results in- 
dicate a non-analytic decay of this mean-field amplitude 
at p — *■ 0, specific for the considered ID system. It is of 
interest to compare this to a power law divergence in p of 
the mean-field amplitudes for systems where a second or- 
der transition at non-zero temperature exists in absence 
of long-range links • 



FINITE SIZE SCALING 
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^ Sum 






• Analytic 
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Figure 6: Comparison between the Monte Carlo simulation 
for a sample of A*' = 8192 and p = 1/4, the numerical sum 
Eq. 11211 and the analytical result (because of the small sizes 
it was necessary to include corrections of 0(A~^/*) to the 
finite size scaling). 



defines the contributions ^ 0{N ^ hiN) into free energy 
/, with 



g{l,u) 



2(1-0 



(2-/(2 + u))(2-u) 



(36) 



and 



It was shown that the self-averaging property invoked 
for our calculation of ^{L, M) is true in the thermo- 
dynamic limit, and fiatness of the histograms in the 
multicanonical simulations suggests that the calculated 
ri(L,M) is accurate (see Fig. (|SJ). Nevertheless, there 
are visible deviations from fiatness near the edges of the 
spectrum, which diminish with growing system size A''. 
One can therefore ask whether our solution also contains 
the correct finite size scaling properties of this model. 

To answer this question, the numerical sum of Eq. Ijl2|l 
was performed for different temperatures. The factorials 
were substituted by the StirHng's approximation and the 
specific heat was then obtained by numerical differentia- 
tion. The comparison between the two independent cal- 
culations is presented in Fig. Despite the fact that 
the Monte Carlo simulation was made for a unique sam- 
ple (without disorder averaging), the agreement between 
these calculations is very good. 

To derive the analytic form of the finite size scaHng 
function, we approximated the sum in Eq. (jl2|l by an 
integral in the intensive variables / and n. Then, after 
the change of variables {l,n) {l,u) where u = n/l, 
the partition function for a finite size A^ of the sample 
becomes 



Zm = exp[-/?Ar(/o + fa)]ZFSS 



(34) 



where 



h{l,u) = 2(1-0 In 
- (2-M)nn 



l-l 

1 - /* 

2 - u 



+ (1 - ul)\n 



l-ul 
1 - u*l* 



l{2 + u))\n 



\u* 
2-l{2 + u) 
2-l*{2 + u*) 



(37) 



The function h(l,u) has a maximum at {l*,u*). When 
T T-, we have /* ^ 0+ and u* 2/[l + 
coth(/3c/)/2] < 4/3. 

The integral in the Eq. (I35|l is dominated by the vicin- 
ity of the maximum. To obtain the leading order terms in 
1/A^ and t we may expand h(l,u) around the maximum. 



*\2 



h{i, u) w -ci{i - ry - C2i{u - u*) 



(38) 



with definite constants ci,2 ~ 0(1), and replace g{l,u) 
by 5(0, u*). With suitable change of variables we obtain 



(C) 



ki 



dx 



(2/u*-l)y* 



dye 



y* 



(39) 



were C — tN^/^, y* cx A^^/*, and the constants ki^2&re 
related to ci,2- Since the finite size scaling limit is N —* 
00 and t — > (at fixed C < 0), the limits of integration 
in y tend to infinity and we get 



Zfss = ^ dl dugil,u)exp{pNh{l,u)) {35) Z^ssiC) ^ ^ / 

27rp Jo Jo f^i Jo 



^-{x + k2Cf^X 



(40) 
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This leads to a correction in the additional free energy: 



fa ~* fa~Tc-^ lnZjr55(C), 



(41) 



and, since fa is also proportional to t^, the scaling form 
for the specific heat becomes: 



Cn{T) 
Coo{T) 



(42) 



Thus the specific heat curves for finite size systems, 
Cn(T), when scaled by Coo(7')and plotted as a function 
of tN^/^, should collapse to a single curve. The results 
of the Monte Carlo simulations are consistent with this 
prediction (see the inset of Fig. ©). An excellent accor- 
dance between the analytic behavior and the results of 
direct summation in Eq. Ijl2ll and of MC simulations is 
shown in Fig. El 

A similar scaHng is observed in the susceptibility (not 
shown here) and has been observed by other authors in 
ID and also in 2D and 3D 0, (where iV L'^, is 
the number of spins, not the linear dimension L of the 
lattice). 

This is the expected form of scaling for a situation in 
which the dimensionality is greater than the upper crit- 
ical dimension and hyperscahng is violated [23l- It is 
observed in all these small world models for any dimen- 
sion of the underlying regular lattice [l^ IT^. A similar 
steepest descent solution probably appHes also in all these 
cases. 



VI. RELATION TO BETHE LATTICE 
APPROACH 

The local environment of a spin in our model looks 
like the Cayley tree in Fig. ijjjl. The vertical Hnks are 
shortcuts (of strength /) and the longer ones segments 
of the ID chain, containing d = l/2p links of strength J. 
Since a shortcut from a given spin has an equal chance of 
linking it to anywhere in the lattice we do not expect to 
find closed loops until we go OjlnN) links away. Based 
on this insight, Dorogovtsev et. al developed a de- 
scription of the Ising model on such lattices based on the 
Bethe lattice solution. Their detailed results (namely for 
Tc) are not directly applicable to our lattices (which are 
not maximally random because of the strong correlation 
between shortcut sites). But their insight certainly is, 
and, since our results are not based on the Bethe lat- 
tice solution, and our lattice has a well defined thermo- 
dynamic limit, it is interesting to consider the relation 
between the two approaches. 

In a Bethe lattice there is only one path to Hnk two 
spins, as a result the correlation function is of ID char- 
acter. Therefore, the correlation between any two spins 
(cTofr) decays exponentially with the distance at any fi- 
nite temperature (even at T^). But the number of spins 
Or at a given distance, r, from a given one, Ns{r), grows 




Figure 7: The Bethe lattice that describes the local envi- 
ronment of any spin in the model: shortcuts (double lines) 
between the sites divide the Ising chain into equal segments 
with d (here d = 3) chain bonds (single lines). 



exponentially with distance, not as a power law, as in 
a regular lattice. The function (aocrr)Ns{r) has a decay 
length that diverges at the ordering temperature of the 
Bethe lattice. It is straightforward to derive the exact 
transition temperature of a Bethe lattice of coordination 
q, PcJ — In {ql {q — 2)) /2, from this condition. 

In our model the chemical distance between two spins 
can be taken as L' -I- M where M is the number of short- 
cuts and L' > M the number of ID chain segments with 
d bonds each, which connect these spins. Then the ID 
correlator is 



(o'oO'(L',M)) = e 



L'd\ntj+M' \nti 



(43) 



The total number of spins with this separation from a 
given one, is 

7V,,(L',Af) = 2*^Cf;, (44) 

therefore 

(aoa(i, M))iV.(i', A/) = e-«(-'^)(^'+^^), (45) 

where the decay constant k{x,T), with x = L'/{L' + M), 
is 

k{x,T) = (l-a;)ln(l-.T) + (2x-l)ln(2a;-l)- 
(1 — x) In 2 — x In a; 

-xh-if'j - {1 - x)lnti (46) 

Minimizing ^(a;, T) with respect to x, we find that this 
minimum value k{T) = minx k(x, T) decreases with tem- 
perature and turns zero just at T = Tc given by Ea. lj30ll . 
At any temperature above Tc the function in Eq. El de- 
cays exponentially and there is no possibility of long 
range order. In the language of the Bethe lattice, the 
occurrence of an extra term in the free energy below Tc, 
expresses the effect of boundaries, which is never negligi- 
ble, no matter what the lattice size is, when the function 
in Eq. 03 does not decay with distance. 
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The finite size scaling for this Bethe lattice again cor- 
responds to the above referred situation when the hyper- 
scaling relations are violated and the finite size correc- 
tions are not determined by the length that characterizes 
the decay of correlation functions. In fact, we have ar- 
gued above that the spin-spin correlation function decays 
exponentially with the distance measured on the equiva- 
lent Bethe lattice, i.e. with the chemical distance on the 
original lattice. Thus, for the instance of / = J, we have 
^ = — 1/ Inij and 



(47) 



It is well known that the chemical distance between any 
two randomly chosen spins is of order C(lniV) so the 
correlation between two spins at a distance of order TV 
along the ID chain is at least of order N~^^^] it does 
not decay exponentially with N as the lattice and the 
distance between spins grow. 

One interesting question that remains unanswered is 
whether one can modify the model in order to effectively 
be at or below the upper critical dimension, and therefore 
observe a non mean-field behavior. 



VII. CONCLUSIONS 

In summary, we have been able to derive an exact so- 
lution of an Ising model on a lattice with long range dis- 
ordered interactions. This solution expresses the free en- 
ergy in terms of the density of states as a function of two 
macroscopic variables of order 0{N), which therefore is 
self-averaging. Hence no disorder averaging is required in 
this approach. We obtained the thermodynamics in the 
N ^ oo limit and also the finite size scaling behaviour. 



VIII. APPENDIX 

The model in which the positions of the shortcut sites 
are randomly chosen can be solved along the same lines 
of the model considered in the main text. We denote the 
coordinates of the shortcut sites, in increasing order, by 
{(ji, 92, • ■ • , q2pN} with qi > (allowing for q2pN = 0) . 
The distances between consecutive sites are 



qt-qi-i, i^l, 
N ~ q2pN + qi 



(48) 



so that Y^'^=i di = N . If we choose M bonds (those for 
which 6=1) from the total of pN , the corresponding 
shortcut sites {ri, r2, . . . , r2M} C {gi, 92, • ■ • , q2pN} will 
divide the lattice into 2M segments of lengths h,. . . , hu 
where 



k = Ti- ri_i, i ^ 1, 
h = N - T2M + ri. 



and YA=x h = ^. Then f)(L, M) is the number of possi- 
ble choices of the M bonds such that 



We define 



^2 + ^4 + ^6 + • • • + ^2M — 



n{L,M) = n{M)P{L\M), 



(50) 



(51) 



where Q.{M) is the number of choices of M shortcuts with 
b = 1 from a total of pN {n{M) = C^f ), and P{L\M) 
is the probability that any such choice of M bonds will 
select L chain bonds with 6=1. 

It should be stressed again that, for a given reaHza- 
tion of disorder, this probability must be calculated in 
the event space consisting of the choices of M shortcuts 
from the specific set of pN random shortcuts. However, 
it follows from the self-averaging property in the thermo- 
dynamic limit that any statistically significant configura- 
tion of shortcuts leads to the same probability P{L\M). 
In that case we can calculate it, enlarging the space of 
events to include all the configurations of shortcuts. 

We are therefore led to ask in how many ways one can 
choose 2M sites, {qi, q2, . . . , (72A/} C {1, . . . ,N — 1}, such 
that the sum of even lengths in this series is ^2 + '4 + ^6 + 
... -I- I2M = I- We have seen above that it is given by 



^N-L^L-l 



(52) 



Since the total number of such choices is C^j , we have 



P{l\s) — C^j ^ C^.j\/ C2M^ 



so that 



n(L,Af) = C?f<r^C,Y-i/C: 



2M- 



(53) 



(54) 



With the same definitions as above we get 

s{l,n) = — pn(lnri/4) — j)(l — n) ln(l 

+ (1 - 2pn)ln(l - 2]jn) 

— (1 — / — pn) ln(l — / — pn) 

— (Z — pn) ln(Z — pn) 

+ (1 - 01n(l - + 'ln^ 

The equation for the transition temperature is 

tjil + Apti) = 1 

with the characteristic limits 

2J . , J 



n) 



(55) 



(56) 



Tc = 
T, = 



ln(l/2p)' 
2J 

In {Tj2pl) 



p < 1, I > 



pi < J. 



ln(l/2p) ' 



(57) 



(49) 



The resulting thermodynamic behavior is essentially the 
same as in the model considered in the main text. 
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